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Summary 

-^This  final  report  describes  the  completed  research  of 
the  coffee  policy  simulation  project.  It  briefly  reviews 
progress  leading  to  the  policy  simulation  stage,  describes 
the  policy  simulation  procedure  itself  and  presents  some 
simulation  results.  The  report  traces  research  activity 
that  has  taken  place  between  August  12  and  September  14,  1981, 
and  consists  of  the  following  parts:  model  estimation  and 
testing,  model  base  forecasts,  and  policy  simulation  of 
quota  levels.  The  theoretical  specification  adopted  for 
the  coffee  model  that  provides  the  basis  for  the  simulations 
is  given  in  the  initial  project  report:  j*Policy  Simulation 
of  the  International  Coffee  Economy:  Model  Description^^t 
Model  Estimation  and  Testing  \ 

The  estimated  or  empirical  model  that  has  been  derived 
from  the  theoretical  model  specification  is  given  in  the  second 
project  report:  "Policy  Simulation  of  the  International 
Coffee  Economy:  Interim  Report."  Missing  from  that  report 
are  the  final  variable  adjustments  to  individual  equations. 
These  are  listed  below  incon^lete  of  the  underlying  statistics 


^Visiting  Professor,  Massachusetts  Institute  of  Technology. 

This  research  is  being  funded  under  a  Department  of  State  PY-81 
external  research  contract. 

DTIC  ' 

8S01  11  088 

H 


Of  parameter  significance,  goodness  of  fit,  and  auto¬ 
correlation.  The  list  of  variables  employed  in  the  model 
is  given  below  in  Table  1,  and  equation  statistics  may  be 
found  in  Appendix  1. 


EXPORTS  BY  COUNTRY  AND  REGION 


brazil: 

exbraeq:  equation 

1>EXBRA=QB^C0NBRA-C  I EBRA- 1  EBRA\  1 ) 

l>AHBR2i  +  <0.5§9341>«AHBRA\1  +  <20.6915>»‘PR1\7  ii 

2>p+^<J^^74394>’<CAHBRA\1-C<281.225>+<0.559341>«AHBRA\2+<20.6913>>‘  SS 

HAS  BEEN  CORRECTED  FOR  FIRST-ORDER  AUTOCORRELATION) 

QBRAEQ:  EQUATION 
r>QBRArAHSM^YLOBRA 

iebraeq:  equation 

5>+<4490 . 96>»DIBRA\1))  ^ 

UEBR>^Q  HAS  BEEN  CORRECTED  FOR  FIRST-ORDER  AUTOCORRELATION) 

YLDBRA,  CONBRA  =  EXOGENOUS 

COLOMBIA: 

excoleq:  equation  , 

l>EXCOi:rQCOL-CONCOL-C  IEC0L-IEC0L\1) 

?>^o8i  <3i237?55>  +  <1.03342>’«QOOL\1  +  <11.14|6>XPR1\1  GG 
CQC^Q  HAS  BEEN  CORRECTED  FOR  FIRST-ORI^  AUTOCORRELATION) 


FIRST- 


CONCOL,  lECOL  s  EXOGENOUS 
OTHER  SOUTH  AMERICA: 

EXS/WTHEREQ:  EQUATION 

1>EXS/V10THE^QS^1>ER-CONS/W1>CR-C  IESAM01TER-  lES/WTHERX  1) 
QS/^THEREQ:  EQUATION 

r>QSAMQTHEfc  <M0.026>  <0.762773>>*QSAMOTHER\l  +  <5.70501>’«PR1\1  GG 

2>>  <388.521>*DS^T«R  - 

C0N5AM0THER,IES/^TICR  s  EXOGENOUS 
AFRICA: 


>»  GG 


f^!l^(^R-ClEAFR.IEAFRM) 
CONRPR,  lEAPR  s  EXOGENOUS 


NORTH  AMERICA: 

EXN/^Q:  EQUATION  ,  ,  ^ 

1>EXN^QN^C0NNAM-C  I ENAM- 1  ENAM\  1) 

?>QNA&  <?112l55>  +  <0.708054>’«QNAM\1  +  <20.5825>’‘PR1\1  M 
2>  +  <1185.16>’«DNAM 

CONNAM^ lENAM  =  EXOGENOUS 

ASIA  AND  OCEANIA: 

EXASIAtOCEEQ:  EQUATION 

1>EXASIA(XE=QA^A80CE-C0NASIAS0CE-(IEASIA80CE-IEASIA80CE\1) 
QASIA&OCEEQ:  EQUATION 

1>QASIA80CE=  <^0.7836>  +  <0.9724S7>’*QASIAtOCE\l  +  <11.2852>«PR1\1  88 
2>  -  <0.392518>’«(QASIA80CE\1-C<-60.78J6>+<0.972497>’«QASIA80CE\2+  88 
3><11.2852>«PR1\2J) 

CQASIA80CEEQ  HAS  BEEN  CORRECTED  FOR  FIRST  ORDER  AUTOCORRELATION) 

C0NASIA80CE,IEASIA80CE  s  EXOGENOUS 

WORLD 


exworldeq:  equation 

l>EXWORL^EXBRAfEXO 


XCOL-«-EXNAM4-EXSAMOTHER-hEXAFR-i€XASIA80CE 


lEPROOTOTEQ:  EQUATION 
1> lEPRODsIEPROOXBRAf  lEBRA 

ieprodxbra  s  exogenous 


IMPORTS  BY  COUNTRY  AND  REGION 


UNITED  STATES: 


MUS3EQ:  EQUATION  ,  „  >  . 

1>MUS=  <28288. 3>  -  <49.6904>>*PR1  -  <3.58319>’*GNP$75US  +  <1.04776>’‘  88 
2>CIEUS-IEUS\1  ) 

EUROPE: 

l>MEuSi  <?OoJ!35>  -  <63.2797>>‘PR1  +  <17.8760>’«GDP$75EUR 
REST  OF  WORLD: 


MOTHEREQ:  EQUATION  .  «  „  . 

l>MOTHEfc  <r66.851>  +  <4.87825>»W$75EUR  +  <1176. a/. 


WORLD: 

MWORLDEQ:  ELATION 
1>MN0RLD=MUS+M^-W 


14MOTHER 


PRICES 


J>?^MWOR^^ljS8LO^J)ilEPROO^?+lxi£RLD)))  +  <18. 

neO'.  EQUATION 

i>prspRp«cptus 

<1.05888>*P1  > 


PCUATEQ:  EQUATI 
l>PCUATs  <07674 
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TABLE  1 


LIST  OF  COFFEE  MODEL  VARIABLES 
(As  utilized  in  the  econometric  specification) 


Symbol 


Identification 


QBRA 

QCOL 

QSAMOTHER 

QAFR 

QASIA60CE 

QNAM 

QWORLD 

EXBRA 

EXCOL 

EXSAMOTHER 

EXAFR 

EXASIA&OCE 

EXMAM 

EXWORLD 

CONBRA 

COM^OL 

CCmSAMOTHER 

COHAFR 

COMASIA&OCE 


Endogenous 

Brazil  coffee  production 

Colombia  coffee  production 

Other  South  America  coffee  production 

Africa  coffee  production 

Asia  and  Oceania  coffee  production 

North  America  coffee  production 

World  coffee  production 

Brazil  coffee  exports 

Colombia  coffee  exports 

Other  South  America  coffee  exports 

Africa  coffee  exports 

Asia  and  Oceania  coffee  exports 

North  America  coffee  exports 

World  coffee  exports 

Brazil  coffee  consumption 

Colombia  coffee  consumption 

Other  South  America  coffee  eonsuuqption 

Africa  coffee  consumption 

Asia  and  Oceania  coffee  consumption 

North  America  coffee  oonsumption 


CONMAM 


PICA76 


PBRICA 


PCOL 


PMIA 

PRl 

PI 

P6UAT 

NOS 

MEOR 

MOTHER 

MNORLD 

lEBRA 

lECOL 

lESAMOTHER 

lEAFR 

IEASIA&OCE 

lEHAM 


ICO  composite  Indicator  price,  1976 
Agreement  (unweighted  average  of 
r^ustas  and  other  mild  arabicas) .  This 
price  series  may  have  to  be  adjusted 
to  reflect  the  new  agreement  formula. 

Unwashed  areUsicas  price  (Brazilian, 
Santos  No  4) .  This  price  series  needs 
to  be  replaced  by  a  better  price  such 
as  U.S.  unit  value  imports. 

ICO  Colombiaui  mild  arabicas  price 
(Colombian  Mams) .  This  price  series 
needs  to  be  replaced  by  a  better  price 
such  as  U.S.  unit  value  imports  or 
Guatonalan  prime  washed. 

ICO  other  mild  arabicas  price  (El 
Salvador,  Central  Standard,  Guatemalan 
Prime  Washed,  Mexico  Price  Washed) . 

U.S.  Unit  Import  Value  (Deflated) 
for  coffee 

U.S.  Unit  Import  Value  for  Coffee 

Guatemala  prime  washed  price 

United  States  net  coffee  imports 

European  gross  coffee  imports 

Rest  of  World  gross  coffee 

World  gross  coffee  imports 

Brazil  coffee  inventory  (end  of  year) 

Colombia  coffee  inventory 

Other  America  coffee  inventories 

Africa  coffee  inventories 

Asia  and  Oceania  coffee  inventories 

North  America  coffee  inventories 

Coffee  inventories  held  by  producers 
other  than  Brasil 


lEPRODXBRA 


lEPROD 

lEUS 

AHBRA 

YLDBRA 

QINV 


Total  producers  coffee  inventories 

United  States  green  coffee  inventories 

Brasil  acreage  harvested  (ha) 

Brazil  coffee  yield  (60  kg  bags/ha) 

Inventories  accumulated 

as  a  result  of  quota  operations 


Exogenous 


GNP$75US 

6DP$75EUR 

T 

j^quota 

CPIUS 

DPR 

DIBRA 

RSTUS 

ACCUS 

DCUS 

DAFR 

DCOL 


GNP  in  United  States  at  constant  market 
prices 

GDP  in  OECD- Europe  \t  constant  market 
prices 

Time  trend  variad>le 

ICA  coffee  export  quota  for  world 
or  for  individual  regions 

United  states  Consumer  Price  Index 

Dummy  variable  for  PRl  based  on  ICA 
quotas  and  1977  member 
disruption 

Dummy  variable  for  extremes  in  Brazil 
coffee  inventories 

Dummy  variable  for  United  States 
reported  green  coffee  roastings 

Dumniy  variable  for  United  States 
apparent  green  coffee  roastings 

Dusnsy  variable  for  extremes  in  United 
States  coffee  roastings 

Dunsqf  variable  for  Africa 
production  cycle 

Duasv  variable  for  Colombia 
production  cycle 


ONAN 


Dmmy  variable  for  North  America 
production  cycle 


7 


DSAMOTHER  Dummy  variable  for  Other  South  America 

production  cycle 

DMI  Dummy  variable  for  extremes  in  rest 

of  world  imports 


The  testing  of  the  model  required  that  the  estimated 
and  actual  values  of  the  endogenous  v^u:iable6  in  the  model 
be  reasonably  close  over  the  saunple  period  of  model  estimation, 
1960-80.  This  closeness  or  accuracy  can  be  measured  in  a 
number  of  ways,  the  most  simple  one  being  the  mean  average 
percent  error  (NAPE) : 


X  ±00% 


where  «  estimated  value  of  a  variable  in  time  period  t, 

A^  »  its  actual  value,  and  n  ■  number  of  time  periods.  Tad>le  2 
shows  the  MAPE  for  different  versions  of  the  coffee  model 
tested.  The  final  model  selected  HISTSIM831  shows 
reasonably  good  accuracy  with  less  than  10  percent  error  for 
most  variables.  Brazilian  exports  proved  as  exception  at 
19'. 0  percent,  largely  because  of  the  model's  Inability  to 
simmulate  erratic  government  intervention  policies. 


Table  2 


MEAN  AVERAGE  PERCENT  ERROR 
FOR  COFFEE  MODEL  ENDOGENOUS  VARIABLES 
OVER  THE  SAMPLE  PERIOD 
1960-1980 


Endogenous 

MAPE 

MAPE 

MAPE 

MAPE 

MAPE 

Variables 

Version 

1  Version  2 

Version 

3  Version 

4  Version 

QBRA 

6.49 

6.46 

8.49 

9.70 

9.81 

QCOL 

3.46 

3.18 

3.62 

4.79 

4.83 

QNAM 

6.94 

6.87 

7.46 

8.13 

7.74 

QSAMOTHER 

4.96 

4.64 

5.56 

6.13 

5.85 

QAFR 

7.91 

7.78 

8.25 

4.53 

4.52 

QASIA&OCE 

4.65 

4.93  ' 

5.95 

7.10 

6.64 

QWORLD 

3.46 

3.12 

2.75 

2.46 

2.39 

MWORLD 

2.45 

2.28 

2.49 

2.85 

2.68 

MUS 

3.99 

3.18 

3.26 

3.50 

3.41 

MEUR 

2.87 

2.73 

2.84 

3.05 

2.98 

MOTHER 

5.57 

5i57 

5.31 

5.31 

5.31 

lEBRA 

12.14 

12.75  * 

31.26 

41.85 

9.12 

lEPROD 

6.72 

6.96 

12.55 

16.30 

4.66 

EXBRA 

14.91 

16.06 

16.71 

18.98 

EXWORLD 

7.40 

7.19 

5.52 

4.80 

4.77 

PRl 

10.30 

12.08‘ 

11.36 

10.76 

11.07 

AHBRA 

ACCUS 

PICA76 

6.49 

6.46 

2.52 

8.49 

9.70 

9.81 

12.41 

1)  Version  WfeblTOO  -  ARl  equations 

s) Version  AU92OA  -  PRl  equation  uses  EXIfORLD  rather 
than  QWORLD 

>) Version  NEWHIST827  -  With  new  USDA  data 
*)VersionHISTSIM828  -  Coffee9HISTl  with  all  new 
equations  required  by  USDA  revisions 
•) Version  HISTSIM831  -  Coffee§RISTl  with  add 
factor  for  lEBRA 


MODEL  BASE  FORECASTS 


Problems  of  Model  Forecasting.  To  evaluate  the 
impact  of  the  coffee  export  quota  levels  over  the  next  five 
years,  a  base  forecast  extending  from  1981-1986  had  to  be 
prepared.  Three  problems  in  particular  had  to  be  solved. 

First,  the  exogenous  variables  in  the  model  such  as 
GNP  and  CPI  had  to  be  forecast  over  the  same  period.  This 
was  accomplished  largely  using  forecasts  in  existing  data 
banks,  mainly  those  of  Data  Resources,  Inc.  Other  exogenous 
variables  to  be  forecast  included  consumption  and  inventories 
in  the  producing  regions. 

Second,  the  special  class  of  exogenous  varieUbles,  the 
dummy  variables,  had  to  be  extrapolated  into  the  future.  This 
was  accomplished  by  careful  analysis  of  foreseable  market 
conditions,  as  explained  in  the  next  section. 

Third,  the  slight  differences  between  FAO  demand  data 
and  USDA  supply  data  had  to  be  reconciled  in  the  forecast 
period  to  facilitate  export  quota  operations  with  the  model. 

This  was  accon^lished  by  equating  imports  with  exports  at  the 
world  level.  A  residual  adjustment  was  then  made  and  allocated 
to  other  importing  regions. 

Model  Adjustment.  The  selection  of  values  for  the 
dummy  variables  in  the  model  followed  the  perception  of  foreseeable 
market  conditions.  These  conditions  in  turn  have  been  checked 
with  expert  opinion  so  that  a  realistic  model  forecast  oould 
be  proposed.  The  variables  together  with  their  values  are 
reported  in  Table  3.  Below  a  rationale  is  presented  for  each. 


Brazxl.  The  unexpected  large  production  of  32,000,000 
bags  for  Brazil  in  1981  required  that  several  related  variables 
be  adjusted.  First,  an  absolute  increment  of  5,000,000  bags 
(QBRA)  was  added  to  the  model's  prediction.  Second,  crop  yields 
that  were  predicted  exogenously  were  increased  from  9  to  11 
bags/hectare  for  1981.  Brazilian  inventories  (lEBRA)  also 
were  decreased  and  placed  in  sales  to  prevent  overaccumulation. 
Finally,  area  harvested  (AHBRA)  as  predicted  by  the  model  grew  too 
sharply  and  the  area  was  decreased  to  mahe  the  prediction 


more  realist!  v. 


TABLE  3 


EXOGENOUS  VARIABLE  ADJUSTMENTS 
FOR  THE  FORECAST  PERIOD,  1981-86 


DAFR  DCOL  OCUS  DIBRA  DMI  Dl^  DPR  DSAMOTHER 


P  8;I8S  8:1118  8; 

1  83  0.7|0  0.000  0. 
1§84  0.750  0.500  0. 
Ills  0.750  0.000  g. 
1986  0.750  0.500  0. 


tftmA 


8QBRA 


8IEBRA 

-’'1 

8 

0 


8:188  .8: 
0.500  0. 
0.500  0. 
o.sog  g. 
0.500  0. 


8:188 

0.500 

81 
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Colombia.  The  dummy  varied>le  reflecting  cyclical 
fluctuations  in  the  coffee  tree  production  cycle  DCOL  was 
changed  from  0-1  variation  to  0  -  0.5  variation  to  reflect 
the  declining  influence  of  this  cycle  in  Colombia. 

Africa y  South  America  and  North  America.  The  cyclical 
production  dummies  were  changed  from  0-1  variation  to  a  constant 
0.5  value  to  reflect  the  declining  influence  of  this  cycler 
i.e.r  DAFRr  DNAMr  and  DSAMOTHER.  The  increase  in  the  value  to 
0.750  in  later  years  for  Africa  reflects  the  cumulative  effect 
of  increased  tree  plantings  on  African  production. 

Europe.  Because  coffee  consumption  in  Europe  is  believed 
to  be  near  saturation  level r  the  model's  prediction  which  did  not 
include  this  factor  were  too  high  and  a  downward  adjustment 
has  been  made.  See  MEUR. 

Prices.  A  dummy  variable  DPR  has  been  used  to  account 
for  direct,  unusual  frost  effects  on  coffee  prices.  An  adjustment 
is  shown  in  1982  to  reflect  the  Brazilian  frost  condition 
carrying  over  from  1981. 

Base  Forecast  Validation.  Given  the  above  model 
adjustments,  a  base  forecast  has  been  produced  which  provides 
the  "most  likely"  coffee  market  scenario  under  free  market 
conditions,  i.e.,  with  no  international  coffee  agreement 
Including  export  quotas  in  effect  during  the  forecast  period. 

The  export  quota  simulations  are  then  tested  with  this  scenario 
as  the  basic  market  outlook.  The  base  forecast  for  the  major 
endogenous  variables  is  summarized  in  Table  4. 

i 


-12- 


Teible  4 

COFFEE  MODEL  BASE  FORECASTS 
1981-1986 


3RA.SET2 
:0L.SET2 
,..-A.SET2 
5SAMOTHER.SET2 

3afr.set2 

^IA80CE.SET2 


1981 

m±. 

18^894.578 

9^522.480 


1982 


1983 


1984 


1985 


14.944 
28 


mm  wp 

mm  mm  ‘m-Jf 


3;255!?? 
20;329.70 


;iil:^25 


1986 


3;20f.9|j 
21,716.186 
10 J 141. 749 


CONBRA 

CONCOL 

CONNAM 

CONSAMOTHER 

CONAFR 

CONASIASOCE 


1981 

8,500.000 

l'850.000 

3'817.000 

1,559.000 

2,697.000 

2,792.000 


1982 

1983 

1984 

8/ 

600.000 

8,700.000 

8,800.000 

i; 

900.000 

1,950.000 

2^000.000 

3/ 

675.000 

3'700.000 

3^725.000 

1, 

575.000 

l'600.000 

1,625.000 

2; 

740.000 

2^790.000 

2,840.000 

2 

800.000 

2,850.000 

2,900.000 

1985 

8,900.000 

2^050.000 

3'750.000 

1,650.000 

2,990.000 

2,950.000 


1986 

000.000 
,100.000 
3;775.000 
1,675.000 
3,040.000 
3,000.000 


I; 


lEBRA 

lEPRODXBRA 

lEPROD 


1981 

19,482.696 

28;431.00Q 

47^913.696 


1982 


1983 


14,362.949 

28^431.000 

42^793.949 


17,415.545 

28;W{.000 


1984 

19,300.200 

28^431.000 

47^731.200 


1985 

21,792.070 

28;431.000 

50;223.070 


1986 

mm 

54;624.331 


EXBRA 

EXCOL 

EXNAM 

EXSAMOTMER 

EXAFR 

EXASIA80CE 


1981 


1982 


1983 


10^120.217 

2,840.070 

15'695.129 

6'827.942 


1984 


14,018.087 

l'803.925 

15,952.002 

6,808.004 


9^150.197 

61913.718 


1985 

'iff 

1,605.667 

17'339.706 

7^033.925 


1986 


mM 

xi'Miii 

7^141.749 


QMORLD 

DCMORLD 

MWORLD 

MUS 

MEUR 

MOTHER 


1981 


1982 


1983 


1984 


Sill 


1985 


ilJ:  1 


1986 


PRl 

PI  , 

PICA76 


1981  1982 


1983 


1984 


.514  4^722  44.38! 

.005  136.031  151.922 

.710  143.270  160.096 


1985 


1986 


itm 

170.621  177.801 


A  typical  problem  of  forecasting  with  a  model  is  that 
the  accuracy  of  the  forecasts  in  a  future  period  cannot  be 
determined,  since  actual  values  are  not  available  for  comparison 
with  forecast  values.  One  approach  is  to  start  the  forecasts 
earlier,  to  "save"  several  periods  for  comparision.  The  approach 
deemed  most  useful  here  is  to  compare  the  model's  forecasts 
with  some  alternative  forecasts  that  could  be  said  to  reflect 
"expert"  judgment.  Below  such  a  comparison  is  made  for  a 
selected  set  of  variables.  The  source  of  the  alternative 
forecasts  are  vinofficial  "expert"  forecasts  by  economists 
from  the  World  Bank  and  other  institutions. 

The  pattern  of  world  coffee  trade  forecast  by  the  model 
is  compared  to  the  expert  forecasts  in  Table  5.  The  expert 
production  forecasts  reflect  a  growth  rate  of  1.4  percent  from 
the  actual  1980  value.  The  model  forecast  of  90,776,000  bags 
for  1986  compares  favorably  to  the  export  forecast  of 
89,943,000  bags.  The  model  forecast  in  the  interim  years, 
however,  reflects  changing  production  conditions.  For  exeunple, 
the  large  1981  and  small  1982  production  forecasts  reflect  the 
Brazilian  influence,  a  bximper  crop  followed  by  a  frost- 
induced  decline.  Other  producers  catch  up  in  1983  but 
relatively  '.ower  prices  in  previous  years  stall  any  further 
production  growth  until  1985  and  1986. 

The  model's  export  forecast  of  62,784,000  bags 
approximates  that  of  the  experts  at  63,752,000  tons  by  1986. 

The  latter  reflects  a  growth  rate  of  1.2  percent  from  the 
1980  actual  value.  Fluctuations  in  intervening  years  reflect 
the  carryover  of  production  from  1981  to  1982  as  well  as 
the  other  stated  production  conditions. 


L980 

L981 

L982 

L983 

L984 

L985 

L986 


PRODUCTI(»«  * 

EXPORTS  * 

Expert 

Model 

Expert 

Model 

Forecast 

Forecast 

Forecast 

Forecast 

82,745 

78,778 

59,850 

55,186 

83,903 

94,850 

60,060 

6.0,651 

85,078 

78,531 

60,783 

62,141 

86,269 

85,363 

61,512 

59,721 

87,477 

84,921 

62,250 

60,147 

88,702 

86,496 

62,997 

60,713 

89,943 

90,776 

63,752 

62,784 

*000/ba98 
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Underlying  the  above  world  production  and  export 
forecasts  are  those  of  individual  countries  and  regions. 

Forecast  comparisons  for  Brazil,  Colombia,  and  Africa  are 
shown  in  Tables  6,  7,  and  8.  The  model's  prediction 
of  Brazilian  production  conditions  reflects  the  buaqper  crop 
and  subsequent  frost  of  1981.  By  1985  the  model  forecasts 
exceed  the  experts'  forecast;  this  largely  reflects  the 
increasing  number  of  coffee  trees  in  Brazil.  The  experts 
forecast  is  based  on  a  growth  rate  of  3.7  percent  from  1980. 
Exports,  however,  are  expected  to  grow  more  slowly:  the 
experts  employed  a  growth  rate  of  3.1  percent  from  1980.  Both 
the  expert  forecasts  and  model  forecasts  are  very  similar 
for  1984,  1985,  and  1986.  The  model  does  not  forecast  a 
greater  export  level;  even  though  Brazil's  production  is 
increasing,  world  in^orts  zure  expected  to  slow  down,  preventing 
a  higher  export  level.  In  addition,  Brazilian  coffee  inventories, 
determined  endogenously,  are  expected  to  increase  over  the  forecast 
period. 

The  model  prediction  for  Colombia's  production  grows 
at  about  the  same  rate  as  the  experts.  The  latter  is  based 
on  a  1.4  percent  growth  rate  from  the  1980  level.  However, 
an  annual  fluctuation  can  be  perceived  because  of  the  cyclical 
crop  production  pattern.  This  same  pattern  is  reflected  in 
Colombia's  exports.  The  grovrth  rate  of  1.0  percent  suggested 
in  the  experts  forecasts  is  less  then  that  of  the  model.  This 


Increase  is  based  on  the  assumption  of  a  relative  increase 
in  the  demand  for  milds. 

The  model  and  expert  forecasts  for  African  production  and 
exports  are  given  in  Table  8.  The  expert  production  forecasts 
are  based  on  a  growth  rate  of  2.4  percent  from  1980  and  the 
export  forecasts  on  3.1  percent.  In  both  cases,  the  forecasts 
are  similar  by  1986.  The  nodel  forecasts  show  more 
realistically  the  impact  of  crop  fluctuations  on  production 
and  exports. 

Table  6 

MODEL  FORECAST  COMPARISON  FOR 
BRAZIL  PRODUCTION  AND  EXPORTS 
1980-1986 


PRODUCTION* 

EXPORTS* 

Expert 

Forecast 

Model 

Forecast 

Expert 

Forecast 

Model 

Forecast 

20,000 
20,740 
21,507 
22,303 
23 


24,872 


18,069 
33,395 
18,017 
25,771 
24,516 
25,3 


27,746 


12,000 


12,755 

13,151 

13,554 

13,979 

14,413 


17,318 

14,536 

14,018 

13,831 

13,922 
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Table  7 

MODEL  FORECAST  COMPARISON  FOR 
COLOMBIA  PRODUCTION  AND  EXPORTS 
1981-1986 


PRODUCTION* 

EXPORTS* 

Expert 

Forecast 

Model 

Forecast 

Expert 

Forecast 

Model 

Forecast 

1980 

14,000 

13,889 

9,750 

9,638 

1981 

14,196 

14,657 

9,848 

10,657 

1982 

14,395 

15,051 

9,946 

12,121 

1983 

14,596 

14,668 

10,045 

11,718 

1984 

14,801 

15,115 

10,146 

12,115 

1985 

15,008 

14,878 

10,247 

1 

11,828 

1986 

15,218 

15,343 

10,349 

12,242 

Table  8 

MODEL  FORECAST  C(»1PARIS0N  FOR 
AFRICA  PRODUCTION  AND  EXPORTS 
1981-1986 


PRODUCTION* 

EXPORTS* 

Expert 

Model 

Expert 

Model 

Forecast 

Forecast 

Forecast 

Forecast 

1980 

19,171 

19,392 

14,052 

14,795 

1 

1981 

18,950 

18,894 

15,441 

15,324 

1982 

19,405 

18,435 

15,920 

15,695 

1983 

19,870 

18,742 

16,413 

15,952 

1984 

20,347 

19,287 

16,922 

16,447 

1985 

20,836 

20,329 

17,446 

17,340 

1986 

21,336 

21,716 

17,917 

18,676 

*000/ba9S 


No  forecast  con^arisons  appear  for  the  other  producing 
regions  in  the  model,  since  expert  forecasts  could  not  easily 
be  assembled  for  these  regions.  The  declining  values  of 
production  and  exports  for  North  America  reflect  the  expected 
continued  dovmward  trend  in  coffee  production  in  this  region 
due  to  coffee  rust  problems,  the  associated  high  cost  of 
production,  and  increasing  political  instzUbility  in  the  major 
producing  countries. 

The  forecast  comparisons  for  world  imports  is  given  in 
Table  9.  The  expert  forecast  is  based  on  a  growth  rate  of  1.2 
percent  from  an  estimate  of  1980  imports.  The  model  forecast 
for  1986  is  65,647,000  bags  compared  to  63,753,000  bags  for 
the  expert  forecasts.  The  world  import  levels  reached  for 
1981  increase  only  slightly  by  the  end  of  the  forecast  period, 
reflecting  the  model's  assumption  of  relative  saturation  of  coffee 
demand  in  the  United  States  and  in  Europe. 

The  final  comparison  is  that  of  prices.  Model  forecasts 
have  been  prepared  for  the  basic  model  price,  the  U.S.  unit 
inqport  value,  as  well  as  the  Guatemala  Prime  washed  price. 

Except  for  1981  and  1982  the  model's  prices  are  above  the 
expert  forecast  price.  The  model  better  reflects  the  changes 
in  production  and  exports  occuring  in  response  to  the  Brazilian 
situation. 
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Table  9 

MODEL  FORECAST  COMPARISON  FOR 
WORLD  IMPORTS* 


Expert 

Forecast 

Model 

Forecast 

1980 

62,000 

62,372 

1981 

60,060 

63,637 

1982 

60,783 

65,158 

1983 

61,512 

65,092 

1984 

62,250 

64,194 

1985 

62,997 

64,998 

1986 

63,753 

65,647 

Table  10 

MODEL  FORECAST  COMPARISON  FOR  PRICES 
(Cvirrent  valuer  1981-1986) 


Policy  Simulation  of  Quota  Levels 

Model  Simulation  Program.  The  theoretical  specification 

# 

of  the  submodel  that  would  predict  the  impact  of  alternative 
quota  level  and  trigger  price  mechanism  policies  has  been  described^ 
in  the  project  papers  cited  earlier,  "Model  Description”  and 
"Interim  Report."  The  translation  of  that  theory  into  an 
effective  agreement  evaluation  submodel  has  been  accos^lished 
by  constructing  an  overall  model  frame«rork.  This  framework 
described  in  Appendix  2  can  be  operated  interactively  with  the 
model  using  the  DRI  network.  The  model  equations  used  for 
the  quota  simulation  2u:e  contained  in  Appendix  3. 

The  program  featuring  the  quota  simulation  frame%rork 
attempts  to  maintain  coffee  prices  within  the  price  range 
specified  by  the  agreement.  Block  I  of  the  model  shown 
in  Appendix  2  decides  whether  coffee  should  be  placed  into  stocks 
or  removed  from  stocks.  In  the  former  case,  the  program 
advances  to  Block  II  of  the  model  and  stocks  reflecting 
differences  between  export  levels  and  quota  levels  are 
stored.  In  the  case  of  higher  prices  and  the  need  to  place 
stocks  on  the  market,  the  program  advances  to  Block  III  of 
the  model.  Coffee  stocks  are  liquidated  to  help  move  prices 
within  the  specified  range. 

Not  included  in  the  Appendix  2  are  a  set  of  additional 


statements  and  equations  that  determine  changes  in  revenue 
resulting  from  quota  operations  for  the  various  exporting  end 
importing  regions  in  the  model. 
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Selecting  Quota  Allocations.  To  perform  simulation 
analysis  with  the  quota  program,  it  is  necessary  to  establish 
quota  levels  for  testing  and  then  to  allocate  them  among  the 
exporting  countries  belonging  to  the  ICA.  The  quota  levels 
to  be  analyzed  are  those  reflecting  the  policy  position  of 
the  U.S.  Government,  some  55-56  million  bags  at  the  world  level. 

To  allocate  these  global  quotas  among  countries,  either 
the  allocation  can  be  given  or  it  can  be  generated  on  the 
basis  of  past  allocations.  The  latter  approach  has  been  employed 
initially.  The  following  listing  shows  the  distribution  of  the  basic 
ICA  annual  quota  of  57,370,000  bags  for  the  crop  year  1980/1981 
together  with  the  distribution  of  non-quota  exports.  Aggregations 
have  been  performed  such  that  the  regions  reported  conform  to 
those  of  the  coffee  world. 


Region 

Basic  Quota 

Non-Quota 

Total 

Brazil 

14.5 

0 

14.500 

Colombia 

9.7 

0 

9.700 

North  America 

10.5 

0.780 

11.280 

Other  S.  America 

2.2 

0.163 

2.363 

Africa 

13.2 

0.984 

14.148 

Asia  and  Oceania 

5.0 

0.373 

5.373 

Total  (000  bags) 

57.364 

Here  non-quota  exports  of  2,300,000  bags  were  allocated 
according  to  the  percentage  distribution  of  the  basic  quota 


among  these  four  regions.  The  total  export  allocations 
shown  in  the  final  column  after  being  converted  to  percentages 
provide  the  basis  for  the  export  allocations  in  the  qn^ota 
program. 

A  simulation  for  the  period  1982-86  was  performed  under 
the  assuo^tion  that  an  International  Coffee  Agreement  (ICA) 
similar  to  the  one  currently  in  effect  (export  quotas  as 
described  in  the  previous  paragraph  and  a  price  range  of 
$1.15-1. 55/lb.)  is  in  operation  over  this  period.  As  shown 
below,  the  simulation  suggests  that  the  ICA  would  not  be 
successful  in  keeping  prices  within  the  specified  range:  the 
price  is  below  the  floor  in  1982  and  above  the  ceiling  in 
1984-86  (the  price  is  within  the  range  in  1983,  as  it  also 
was  in  the  base  simulation) .  This  appears  to  be  because 
the  stocks  accumulated  in  the  (unsuccessful)  attempt  to  raise 
prices  in  1982  are  exhausted  in  1984,  the  first  year  of 
pressure  on  the  ceiling,  leaving  little  in  the  way  of  stocks 
to  defend  the  ceiling  in  1985-86. 

One  point  that  should  be  made  at  this  time  is  that  a 
comparison  of  projected  "free  market"  export  levels  (Table  4) 
with  the  export  quotas  used  for  this  simulation  shows  that 
Brasil,  North  America,  and  South  America  will  likely  not  be 
able  to  meet  their  quotas  1982-86,  while  Coloinbia,  Africa 
and  Asia  and  Oceania  could  easily  exceed  theirs.  This  means 
that  the  pattern  of  inventory  accumulation  and  liquidation 
simulated  by  the  model  may  not  accurately  reflect  the  actual 
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pattern.  The  ICA  has  provisions  regarding  the  reallocation 
of  export  quota  shortfalls,  and  clearly  it  would  be  desirable 
to  include  an  export  quota  reallocation  feature  in  the  next 
version  of  the  coffee  model. 


PRICE 

(PICA76) 


Year 

1982 

1983 

1984 

1985 


Base  Forecast 
$0.98 
1.43 
1.60 
1.71 

7  -  7fi 


Quota  Forecast 
$1.09 
1.55 
1.67 
1.77 
1 


Equation  Statistics 


AHBRAEQ 

LEAST  SQUARES  WITH  FIRST-ORDER  AUTOCORRELATION  CORRECTION 


ANNUALC1961  TO  1980) 
DEPENDENT  VARIABLE: 

20  OBSERVATIONS 

AHBRA 

COEFFICIENT 

STD.  ERROR 

T-STAT 

INDEPENDENT  VARIABLE 

281.225 

227.5 

1.236 

CONSTANT 

1)  0.559391 

0.1301 

9.299 

AHBRAM 

2)  20.6913 

6.896 

3.000 

PR1\7 

0.479394 

0.2763 

1.717 

RHO 

R-BAR  squared:  0.9599 
DURBIN-WATSON  STATISTIC:  1.9325 
STANDARD  ERROR  OF  THE  REGRESSION: 

157.7  NORMALIZED:  0.05739 

lEBRAEQ 


LEAST*  SQUARES  WITH  FIRST-ORDER  AUTXORRELATION  CORRECTION 


ANNUALC1961  TO  1980) 
DEPENDENT  VARIABLE: 

20  OBSERVATIONS 
lEBRA 

COEFFICIENT 

STD.  ERROR 

T-STAT 

INDEPENDENT  VARIABLE 

-17556.0 

2236 

-7.853 

CONSTANT 

1)  0.883839 

0.03300 

26.78 

(IEBRA\l4QBRA) 

2)  9990.96 

1171 

3.835 

DIBRA 

0.908299 

0.2933 

1.678 

RHO 

R-BAR  SQUARED: 

iVS8B^r 


NORMALIZED:  0.06636 
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LEAST  SQUARES  WITH  FIRST-ORDER  AUTOCORRELATION  CORRECTION 


ANNUALC1961  TO  1980) 
DEPENDENT  VARIABLE: 


20  OBSERVATIONS 
QCOL 


COEFFICIENT 

STD.  ERROR 

T-STAT 

INDEPENDENT  VARIABLE 

-1237.55 

408.7 

-3.028 

CONSTANT 

1) 

1.03342 

0.06784 

15.23 

QCOLM 

2) 

11.1456 

5.426 

2.054 

PRIM 

3) 

1437.51 

165.8 

8.669 

DCOL 

-0.631282 

0.2014 

-3.135 

RHO 

R-BAR  SQUARED:  0.9776 

DURBIN^TSON  STATISTIC:  1.8176 

STANDARD  ERROR  OF  THE  REGRESSION:  297.6  NORMALIZED:  0.03313 

QSAMOTHEREQ 

ORDINARY  LEAST  SQUARES 

ANNUALCmi  TO  1980) 
DEPENDENT  VARIABLE: 

20  OBSERVATIONS 
QSAMOTHER 

COEFFICIENT 

STD.  ERROR 

T-STAT 

INDEPENDENT  VARIABLE 

280.026 

348.5 

0.8034 

CONSTANT 

1) 

0.762773 

0.1267 

6.018 

QSAMOTHERM 

2) 

5.70501 

3.501 

1.629 

PRIM 

3) 

388.521 

122.5 

3.170 

DSAMOTHER 

R-BAR  SQUARED:  0.8263 
DURBIN-WATSON  STATISTIC:  2.0513 
STANDARD  ERROR  OF  THE  REGRESSION: 

0^ 


236.1  NORMALIZED:  0.07479 


ORDINARY  LEAST  SQUARES 


0)  ^^OBSERVATIONS 


COEFFICIENT 

STD.  ERROR 

t-sTat 

INDEPENDENT  VARIABLE 

2112.55 

1235 

1.710 

CONSTANT 

1) 

0.708054 

0.1185 

5.975 

0^1' 

2) 

20.5825 

12.24 

1.682 

PRIM 

3) 

1185.16 

576.4 

3.149 

wm 

2.1949 

STAN>ARD  ERROR  OP  TIC  REGRESSION: 


881.6  NORMALIZB):  0.07107 


Ai?pendix  1  (continued) 
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Appendix  1  (continued) 


MEUREQ 


ORDINARY  LEAST  SQUARES 


ANNUALCiaeO  TO  1979) 
DEPENDENT  VARIABLE*. 

COEFFICIENT 

2008.35 

1)  -63.2797 

2)  17.8760 


20  OBSERVATIONS 
MEUR 


STD.  ERROR 

T-STAT 

854.6 

2.350 

11.85 

-5.341 

0.7094 

25.20 

INDEPENDENT  VARIABLE 
CONSTANT 

0 

PRl 

GDP$75EUR 


R-BAR  SQUARED*.  0.9769  ,  „„„ 
DURB IN-WATSON  STATIgIC;  1.87 « . 
STANDARD  ERROR  OF  THE  REGRESSION. 


834.2 


NORMALIZED:  0.03255 


MOTHEREQ 


ORDINARY_LEAST  SQ^RES 


ANNUALC1960  TO  1979) 
DEPENDENT  VARIABLE: 

COEFFICIENT 

166.851 

1)  4.87825 

2)  1176.67 


20  OBSERVATIONS 
MOTHER 


STD.  ERROR 

T-STAT 

552.0 

0.3023 

0.3761 

12.97 

419.4 

2.806 

504.4 


INDEPENDENT  VARIABLE 
CONSTANT 
GDP$75EUR 
DMI 

NORMALIZED:  0.06927 


R-BAR  SQUARED:  0.9031 
DURBIN-WATSON  STATISTIC*.  2.0761 
STANDARD  ERROR  OF  THE  REGRESSION 


PRIQEQ 


ORDINARY  LEAST  SQUARES 


ANNUALC1962  TO  1980) 
DEPEI«)ENT  VARIABLE: 

19  OBSERVATIONS 

PRl 

COEFFICIENT 

STD.  ERROR 

T-STAT 

INDEPENDENT  VARIABLE 

1)  0.218438 

0.08613 

2.536 

PRIM 

2)  28.5251 

3.914 

7.288 

3)  18.8667 

2.850 

6.620 

DPR 

Appendix  1  (continued) 


PICA76EQ 


ORDINARY  LEAST  SQUARES 


ANNUAL(1960  to  1980) 
DEPENDENT  VARIABLE: 

21  OBSERVATIONS 

PICA7S 

COEFFICIENT 

STD.  ERROR 

T-STAT 

INDEPENDENT  VARIABLE 

-0.770841 

3.565 

-0.2162 

CONSTANT 

1)  1.05888  ,  0.04119 

R-BAR  SQUARED:  0.9706 
DURBIN-WATSON  STATISTIC:  1.7458 

25.71 

PI 

DURBIN-WATSON  STATISTIC:  1.7458 

STANDARD  ERROR  OF  THE  REGRESSION:  9.987  NORMALIZED:  0.1392 


PGUATEQ 


ORDINARY  LEAST  SQUARES 


ANNUALC1960  TO  1979) 
DEPENDENT  VARIABLE: 


1) 


20  OBSERVATIONS 
PGUAT 


COEFFICIENT 

0.674788 

1.14411 


STD.  ERROR 
2.774 
0.03451 

Stic:  1.9698 


t-STAT 

0.2433 

33.16 


R>BAR  SQUARED:  0. 
durbin-vRtson  STAT 
STANDARD  ERROR  OF  THE  REGRESSION:  7.597 


INDEPENDENT  VARIABLE 

CONSTANT 

PI 

NORMALIZED:  0.1035 
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LOGICAL  STRUCTURE  FOR  OPERATION  OF  COFFEE  MODEL  UNDER 
INTERNATIONAL  COFFEE  AGREEMENT  (ICA)  CONTAINING  EXPORT 
QUOTAS  AND  PRODUCER  ACCUMULATION/LIQUIDATION  OP  STOCKS 
TO  KEEP  COFFEE  PRICE  WITHIN  RANGE  SPECIFIED  BY  AGREEMENT 

I.  Basic  Coffee  Model  Solution 

SOLVE  (t)  COFFEEMODEL 

IF  Price  less  than  $1.15  THEN  "Export  Quotas  emd  Stock 
Accumulation" 

IF  Price  greater  than  $1.55  THEN  "Stock  Liquidation" 

IF  $1.15  <  Price  <  $1.55  THEN  SOLVE  (t+1)  COFFEEMODEL 


II.  "Export  Quotas  and  Stock  Accumulation" 

QEXWORLD  «  ICA  global  export  quota  (000  bags)  >  a  series 
QEX@Region  »  ICA  export  quota  assigned  to  each  country  or 
region  (000  bags)  «  a  series 
QQ  »  export  quota  cuts  »  V(. 975, .950, .925, .900) 

EXPORTS §Region  ■  QQ*QEX@Region  and 

QINV@Region  =  stocks  accumulated  through  quota  operation 
«  QEX@Region-EXPORTS@Region 
SOLVE (t)  COFFEEMODEL  for  new  EXPORTS @Region 
IF  $1,115  <  Price  <  $1.55  THEN  SOLVE (t+1)  COFFEEMODEL  ELSE 
try  next  QQ  and  new  EXPORTS @Region 

IP  after  all  QQ  have  been  tried  and  $1.15>  Price  THEN  exit 
with  message  "Quota  cuts  and  stock  accumulation  are 
insufficient  to  move  price  into  ICA  range;  this  occurs  in 
year  t." 


III.  "Stock  Liquidation" 


LQQ  ■  rate  of  stock  liquidation  «  V(0.1  to  1.0,  step  0.1) 
EXPORTS §Region  ■  QEX§Region  +  LQQ*SUM (QINV@Regionl 
SUM(QINV§Region)  ■  SUM(QINV@Region)  -  LQQ*SUM(QINVeRe9ion) 

(this  expression  can  be  negative) 

IP  SUM(QINV@ALL  Regions)  >  0  THEN  SOLVE  (t)  COFFEEMODEL  with 

new  EXPORTS ^Region  ELSE  exit  with  message  "Liquidation  of 
stocks  accumulated  under  quota  operation  is  insufficient 
to  move  price  into  ICA  range;  stocks  reach  0  in  year  t." 

IF  $1.15  <  Price  <  $1.55  then  SOLVE  (t+1)  COFFEEMODEL  ELSE 
try  next  LQQ  (eventually  SUM(QINV$ALL  Regions)  will  reach  0) 


Model  Used  for  Quota  Siamlation 


qcoleq:  equation 
r>QCOi:=.  <-1237.5 
2>  +  <1437.51>^I> 
3>QC0L\2+<1 1.1451 


542>*  U 


qnameq:  equation 

r>QNAfe  <2112.55>  +  <0.708054>»QN^\1  ♦  <20.5825>*PW\1  88 
2>  +  <1185.16>«DNAM 

QASIA80CEEQ:  EQUATION, 

r>QASIA80CE=  <^g.7836>  +  <0.572437>“QASIA80CE\1  ♦  <ll.2852>*pRia  S 


3X11 


80.026>  4  <0.762773>*K}5AH0THER\1  *  <5.70501>»PR1\1  88 
DSAMOTHER 


AHBRAEQ:  EQUATION 
l>AHBRte  <m.225: 


1>AHBRA=  <28l.225>  +  <0. 
2>  +  <0.474394>«CAHBRA\1 
3>PR1\8)) 


225>  +  <p.559341>«AHBRA\l  ♦  <20 .6913>>*PR1\7  88 
>«CAHBRA\1-C<281  .225>+<0 .559541>*AMB^\2-k20T6 


913>*<  88 


qbraeq:  E^TION 

r>QBRAsAHWS«YLDI 


YLOBRA 


?>^rS’<|?!{o7!5?>  +  <0.990399>’*QAFR\1 
2>CCCPR1\8+PR1\9+PR1\10)/3)5  +  <5045.9 


I  +  <31.4871>x  88 
99>»*t)AFR 


QWORLDEQ:  EQUATION 
T>QkK)RLDsQBRAi^Q^-»QNAH«QSAMO' 


ITHER+QAFR4QASIA80CE 


QINVBRAEQ:  EQUATION 

r>OIN^RA=IFTQBRA-CONBRA-ClEBRA-IEBRA\l))  LEQ  QEXBRA  THEN  QINVBRAM 
2>aQBRA-C0NBRA-ClEBRA-IEBRA\l))-QEXBR>O<H)INymRA\l 

?>WNSco8ilF^^ilQCoffa)NCOL-ClE(»L-IEC^  LEQ  QEXCOL  THEN  QINVCOLM 

2>CCVOL-COrCOL-ClECOL-IECOL\l)3-QEXCOL3+QINVCOL\l 

QINVNAHEQ:  EQUATION 

r>^N^%IF  LQNAM-C0(«m-CIEN^IEN4M\1>)  LEQ  QEXNAM  THEN  QINVNAMM 
2>CCQNAM-aMm-ClENm-IENM1\l>)-QE}Q4AH)4QINVNAM\l 

?>mW^ilF^*i®[fSicON^-ClEAFR-IEAFR\l))  LEQ  QEXAFR  T«N  QIWAFRM 
2>CCQAFR-C0NAFR-CIEAFR-IEAFR\1))-QEXAFR)+QINVAFR\1 
QINVASIA80CEEQ:  EQMTION 

r>QiNVASIA80CEsIFX^IA80CE-CONASIA80CE-ClEASIA80CE-IEASIA80CE\l))  LEQ 

3>^^^^L^jNAS?iS^^Si)A8%I^I&lA80CE\l))-QEXASIA80CE>f 
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QIWSAMOTHEREQ:  E(XJATION 
r>QINVS/^THEfcIFXC 


WS/^THEfc  I F  XQSAM0THER-C0NSAM01>CR-C  I  ES/V10THER-IESAMDTHER\  1))  LBQ 

2>QEXSAM0THER  THEN  QINVSAMQTHERM  ELSE  SS 
3>CrQSAMOTHER-Ca 
QINVSAMOTHERM 


^  QINVSAMOl  _ 

iSAMOTHERrCC)NSAMOTHER-ClESAMOTHER-IES/WOTHER\l))-QEXSAMOTHER)+ 


exbraqeq:  equation 

l>E)ffl^IF  C^RA-C0NBRA-CIEBRA-IEBRA\1))  LEQ  (JEXBRA  THEN  EXBRA  ELSE  tt 
2X)EXBRA 


EXCOLQEQ:  EQUATION 

1>EXC^?IF  C^OL-CONCOL-ClECOL-IECOL\l))  LEQ  QEXCOL  THEN  EXCOL  ELSE  U 


EXSAMOTHERKQ:  EQUATION 

1>EXSAM0THER=IF  TOSAMOTHER-CONS/WOTHER-ClES/WTHER-IES^THERM))  LEQ  SS 
2>QEXSAMOTHER  THEN  EXSAMOTHER  ELSE  QEXSAMOTHER 


EXAFRQEQ:  EQUATION 
;5if  CQAFR- 


1>E> _ 

2>QEXAFR 


QAFR-C0NAFR-CIEAFR-IEAFR\1))  LEQ  QEXAFR  THEN  EXAFR  ELSE  U 


EXASIAtOCEQEQ:  EQUATION 

l>EXASIASOasIF  CQ^IA(XE-CONASIA&OCE-ClEASIASOCE-IEASlASOCE\l))  LEQ  U 
2>QEXASIA(OCE  THEN  EXASIASXE  ELSE  QEXASIASOCE 


exworldeq:  equation 

l>EXWORLO=EXBR^XCOL-t€XNAM4€XSAMOTHER<t-EXAFR^XASIASOCE 


EXNAMQEQ:  EQUATION 

1>EXN$€iF  (^NAM-CONNAM-CIENAH-IENAMM))  LEQ  QEXNM1  THEN  EXNAH  ELSE  U 
2>QEXNAM 


<28.5251>’(  SS 

2>CCCMWORLD\l-fMWORLD\25/CIEPROD\i<(€XMO 


IEPROD\l<(€XMORLD)»  *  <18.8667>>(DPR 


iVSSSi  .  3>-<4?  .6904>*PRl-<3 . 58319>’‘GNP$75U&Kl .  04776>>»W 

2>CIEUS-IEUS\1)3  LEQ  Al^EXMORLD  THEN  MUS  ELSE  aI^EXWORLD 


l>MEU^iF  C<2008.35>-<63.2797>’*PR1-K17.876>*GDP$75EUR)  LEQ  88 
2>A2KEXWORLD  THEN  MEUR  ELSE  A2KEXW0RLD 


l>tSoT®SiF^?aJJ?851>+<4.87825>’«GOP$75EUR+<1176.67>>»DMO  LEQ  88 
2>A3’‘EXWORLD  THEN  MOTHER  ELSE  A3>(EXWORLD 


lEPROOTOTEQ:  EQUATION 
1>  I  EPRODs  I  Emoroe^iEBRA 

llama  .  <1.05888>xpi 

MMORLOEQ:  EQUATION 
1>MWORLOsMUS«MEUR4MOTHER. 

{>8!S?s  -  <40.0123>>4»R1  -  <2.89770>»«CNP$75US  88 

2>  *  <987.383>^DCUS 


TION 


INVCOL4QII^NAMfQIfNSAMOTHEI4fQI»MAPI8>QI»^ASIASOCE 


impend ix  3  (continued) 


!sm  iOTi>  ♦  <1.35736>»«P1 

pcoleq: 


i>PcoS 


JION 

0259>  +  <1.13986>*P1 


<1.13934>»‘P1 


POTHERMILDEQ;  EQUATION,,  „ 
1>P0THERMIL0=  +  <1.09966>*P1 

t5m!iiOTSS4>.<l.UU2>xpX 


Routines  (Based  on  Appendix  2) 
Used  for  Quota  Simulation 


SOLVE:  ROUTINE 
0.5>DO  FINISHER 

1  >SET  INT=82  TO  86 

2  >ORIGMODELs(X)FFEE8QUOTA2 

3  >D0  START 

START:  ROUTINE 
9*5  I  dated  0INT 
0.75>OISPLAY  "SOLVING 


jTRINGCYEARCSTARTDATECeiNT))),/ 

1  >SOLVE<WARNINGSsNULL>  ORIGMOD^ 

2  >IF  PICA76  LEO  115  THEN  DO  LOW975 

>  ELSE  DO  FINI 


> 

> 

> 

>END 


lER 


THEN  DO  SUMINVSOl 


DECIDE2:  ROUTINE 

1>IF  PIC/^  LEO  115  THEN  DO  LOW975 
2>  ELSE  DO  FINISHER 

ROUTINE 

iTA2=QRIGMgOEL  ^ 
iTSsO^RC^FEE8QUOTA2) 


lALL)  TICN  DO  NEMQINV 


9>rF  sreciFi^iuMB, 

m  RSISIr 
8:} 

ITA2 

5_THEN  do  LOW950 


>SOL^COFf^l 

>0ECI0E 
> 


THEN  DO  SUMII^Ol 
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„  .,30:  ROUTINE 
1>50  FINISHER 
2>L00P_I  BY  QUOTAS  BEGIN 
3^  Isl  #95 
4>  END  „ 

5>SOLVE  COFFEEeQUOTA2 
5>jP  IRACKtowrycm 

»> 


TRACKLOUteVCTRAOCLOW, 

PICATC  LEO  115  THBJ'DO  LOW925 
ELSE  GO  TO  DECIDEl 
9>DECIDE1:IF  PICA76  GEQ  155  THEN  DO  SUMINVSOl 
10>  ELSE  DO  FINISHER 


LOW925:  ROUTINE 
1>D0  FINISHER 
2>LOOP  I  BY,QUOTAS  BEGIN 
!>  1=1^.925 
4>  END  ^ 

5>SOLVE  COFFEEeqiJOTA2 
6>  TRACW.O»iteVCTRAaCLOW..925) 

7>IF  PICS^LEQ  115  THEN  DO  LOW900 
8>  ELSE  W  TO  DECIDEl 
9>DECIDE1:IF  PICA76  GEQ  155  THEN  DO  SUMINYSOl 
10>  ELSE  DO  FINISHER 

LOW900:  ROUTINE 
1>D0  FINISHER 

2>  LOOP  I  BY  QUOTAS  BEGIN 
!>  Irl*.9  ^ 

END 

5>SOLVE  COFFEE0QUOTA2 
|>  TRACkLOWs^TRACKLOW..93 

7>IF  PICA76  LEQ  115  THEN  U 

b^ISPLAY  //•QOoTA  cuts  ARE  INSUFFICIENT  TO  MOVE  PRICE  INTO  ICA  RANGE", 
^j>a.SE  GO  TO  DECIDEl 

T«,  DO  SUMIWSOl 

DOSUMS:  ROUTINE 

EGIN 


4>END 

newqinv:  routine 

UQfNVB^SUMfiBRA 


jS/VOTHER 

5>QlNVAFRsSUM|AFR  _ 
6>QINVASIAG0C&SUM9AS1AG0CE 

SUMINVSOl:  ROUTINE 

I 


11>00  FlNlShCR 


lALL  GEQ  0  THEN 
SPLAY  7,"QU0TA 


.  ^^uSBNM^SUMeS/V10THER-l>SUMBAPR4>SUM8ASIAfiOCE3 
DO  Hicifioi  _  , 

STOCKS  REACH  0",/ 
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ROUTINE 


HIGHOl:  ROUTINE 

1  >COFFEE§(^TA2=COFFEE|gUOTA2  EXCEPT 

2  >EXSAMOTHER()EQ,EXAFRQE^EXASIA(OCEQ 

8  >TRACICHIGH=W.D 

9  >IF  PICA76  GE6  155  THEN  DO  SUMINVSO 

10  >  ELSE  TO  DECIDE2 


>IF  PICA76  GE6  155  THEN  DO  SUMINVS03 
>  ELSE  TO  DECIDE2 


SUMINVSOS:  ROirriNE 
l>SUl«BRA=SU^BRAri .  7 
2>SUM  iC0L=SU^«C0L^ .  7 
3>SUM  iNAMsSU^NAMPt .  7 
‘f>SUM  iSAMOTHEfcSUM@SAMOTHER’« .  7 
5>SUM  |aFR=SUM@AFR’‘,7 
|>SU^«ASIAfiQCE=SUM@ASIj^OCe‘ .7^ 

9>ELSE  DISPLAY  7/'()UOTA  STOCKS  REACH  O", / 

10>DO  FINISHER  ' 

HIGH03:  ROUTINE 

l>COFFEE|qWTA2=COFFEE0QUOTA2  EXCEPT  NLCEXBRAE' 

|>COFFEEeQUOTA2=ORDERCCOFFEE8$JOTA2) 

7>SOLVE  COFFE||qUOTA2 
8>TRACKHIGH=VCTRACKHIGH. .3) 

J>IF  PICA76  GEd  155  TH^  DO  SUMINVS05 
10>  ELSE  TO  OECIOE2 


I0IA2^EXCEPT  NLCEXBRAE 

Bfii^rNiM 


QH03, 


;05:  ROUTINE 
BRA=SUMBBRA?<.5 
COUSUmCOL^.S 
NAMsSUMNArf<.5 
SAMO^RsSUJJSSAMOTHER* .  5 

l/^^i&^BRAfSUMScSLlsU^^ 

UM^L  g1(}  0  THEN  DO  HIGH05 
dTsPLAY77"(5UOTA  stocks  reach  0'% 


SUMINVSOS:  ROUT 
l>SiM  BRAsSlMB 

3>l!»f  SSi^iliSi 
4>Sg«  SAMOThg 

|>iS«ASIAi^ 
7>sumIall^sum 
8>if  sumSall  q 

9>ELSE  DISPLAY 
10>DO  FINISHER 


^^COL+SyMgNAMfSUM0SAMOTHER+SUM8AFR+SUM8ASIAWCE) 
(5UOTA  STOCKS  REACH  0'%/ 


HIGH05:  e 
1>C0FFEE8 

5>COFFEEi 

fef&aMI 

|>IF  Pig 


sCOFFEEl 


EXNAMEi 

ZsORDERC 

IQU9TA2. 


EXCEPT  NLCEXBRAEQH03.EXCQLE(3H03,GG 
^g§ftSM«EQ^«^EXAfRE<lA05,EXASlAite 


iacHi^.52 

^LmOTHERK.3 


SUMI^^S07 


L4>SUMeSM«)THER<fSUM8AFR4>SUNeASIAtOCE) 
EACH  0",/ 
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MIGH07:  ROUTINE 

l>C0FFE^qU0TA2sC0FFEEggU0TA2  EXCEPT  NLCEXBRAEI 

2>EXNAMEQH05.EXSAMOTHE|E()W5.EXAFRE9H05. - 

J>COFFEE0QUOTA2=COFFEE8(^TA2 .  CONCA^ 

i>COFFEE0QUOTA2io^^C6FFeE0QiX)TA2) 

7>S0LVE  CO^Ee2^TA2  ^ 
8>TRACXHIGH=VCTRAaCHIGH. .7) 

9>IF  PICA76  GEO  155  THEN  DO  SUMINVS09 
10>  ELSE^  OEC1DE2 


^EQH05,G( 

1EQH07, 


09:  ROi^INE 
BRA=Sli«BR/^.001 
!COLsSU^iCOL^.001 
NAMnSUWNAMpi.OOl 
SAMOTHEteSU^SAMOTHER>< .  00 1 
AFR=SUM0AFR»«.OO1  _ 

ASI  A«OCE=SUM@ASIMOCEX .  00 1 


J>IF 

)>ELSE  PRINT  SI 
19> 

16>00  FINISHER 


AfiOCE=SUM@ASI>^OCE»‘ .001 

,^SUMgBRAfSUM0COL-i-SUM|NAHfSUM0SAMOTHER-('SUNeAFR-»-SUMSASIAGOCE) 
[all  geo  100  THEN  DO  HIGHOS 
lINT  SUr^ALL 


HIGH09:  RC 
l>COfrEC 
2>EXNAME 


JTINE 

3TA2=C0FFE1 
)7.EXSAMOTHEI 


3>COFFEE0QUOTA2=COFFEE0QUOTA2  CONCAT  _ _ 

6>COFFEE0QyOTA2rORDERCCOFFEEl$JOTA2) 

7>SOLVE  COFFEE|qUOTA2  ^ 

8>TRACXHIGH=VCT^CXHIGH,  .9) 

9>IF  PICA76  GEO  155  TH^  (G 

10>DISPLAY  /."(^TA  STOCKS  EXHAUSTED  BUT  PRICE  STILL  NOT  IN  ICA  RANGE*',/ 
11>PRINT  SUMfiAa 
12>DO  FINIShCR 
19> 

?  PR I NT  EXBRAEQHO 1 , EXBRAECJHO 3 , EXBRAEQHO  5 , EXBRAEQHO 7 , EXBRAEQHOS 
EXBRAEQHOl EQUATION 

^^^C^RA>CIEBRA-IEBRA\1»  LEQ  QEXBRA  THEN  EXBRA  ELSE  U 
^RA-CIEBRA-IEBRAM))  LEQ  QEXBRA  THEN  EXBRA  ELSE  GS 

EXBRAEQH05 EQUATION 

1>^^^  ^^A|CONJRA-C«EBRA-IEBRA\l))  LEQ  QEXBRA  THEN  EXBRA  ELSE  GG 

>CIEBRA-IEBRA\1))  LEQ  QEXBRA  THEN  EXBRA  ELSE  GG 
CIEBRA-IEBRAM))  LEQ  QEXBIM  THEN  EXBRA  ELSE  GG 


